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Abstract 

The intersection numbers of p-spin curves are computed through corre- 
lation functions of Gaussian ensembles of random matrices in an external 
matrix source. The independence of intersection numbers is determined as 
polynomial in p ; the large p behavior is also considered. The analytic contin- 
uation of intersection numbers to negative values of p is discussed in relation 
to SL(2,R)/U(1) black hole sigma model. 
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1 Introduction 



The intersection numbers for p-spin curves appear in the generahzed Kont- 
scvich matrix model [1, 2, 3]. The generating function for p-spin intersection 
number obeys the p-th KdV equation or Gehand-Dikii equation. In a random 
matrix theory, the correlation functions at the edges of the spectrum, where 
one can tune a degeneracy of order p, are expressed through intersection 
numbers [4, 5, 6, 7]. In conformal field theory, the p-spin curve intersection 
theory is related to J\f=2 super conformal minimal theory for Lie algebra Ap_i 
type. It has been pointed out that it corresponds to a gauged Wess-Zumino- 
Witten (WZW) model of SU{2)k/U{l), where /c = p - 2 is the level of the 
Kac-Moody algebra of Lie group SU{2) [8], and is related to SL(2,R)/U(1) 
black hole sigma model when k becomes negative [9] . 

The free energy for the p-spin curve satisfies interesting universal equa- 
tions, such as a string equation, dilaton equation, and WDVV equation, so 
called tautological equations or universal equations, and has been studied in 
the connection to Gromov-Witten theory [10, 11, 12]. Although the p-spin 
curve intersection numbers can be obtained through tautological equations 
in a recursive way, the actual computation for higher genuses is limitted 
[13, 14, 15]. 

In previous papers, we have derived explicit integral formula for the p- 
spin curve intersection numbers of the moduli space „ valid for all order 
of genus g. We have shown that they are obtained analytically for a fixed 
number n of marked points. Our formulation starts from simple Gaussian 
matrix models with an external matrix source and based upon a duality 
relation [5, 6, 7], from which one recovers a generalized Kontsevich matrix 
model. 

The intersection numbers for the spin moduli spaces with n-marked points 
are obtained from the n n-point correlation functions f/(si, Sn) of Gaussian 
random matrices in a scaling limit at critical edges [16, 17]. In a previous 
article [18], we have computed explicitly the intersection numbers of moduli 
space of p-spin curves with one marked point, for arbitrary values of p, as 
polynomials in p. This allowed us to consider continuations in p ; in particular 
the limit p — >■ — 1 exhibits an interesting relation between the intersection 
numbers, and the orbifold Euler characteristics x{Mg,i) — C(l ~ 2gf), where 
(^{x) is the Ricmann zeta function) [20, 21]. 

In this paper, we extend the evaluation of the intersection numbers be- 
yond the one marked point {n — 1) for arbitrary p. The obtained intersection 
numbers are consistent with previously known results [13, 14, 15] for small 
values of p. We pursue the large p behavior, p — )■ cxo limit. The p-spin curve 
intersection theory is equivalent to gauged WZW model. For this gauged 
WZW theory, in which k — p — 2 appears as overall factor, the large k hmit 
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may give a semi-classical solution [9, 23]. In the negative k, the gauged 
WZW model on SU{2)/U{1) is changed to WZW model on non-compact 
SL(2, i?) /[/(!), which is relevant to a black hole a model [9]. We will discuss 
the relation between the intersection numbers and the density of state of 
SL{2,R)/U{1) black hole sigma model [24, 25, 27]. 



2 Generating function for p-spin Intersection 
numbers 

The mathematical definition of the intersection numbers of the moduli space 
of p-spin curves with s-marked points is given by[3] 

< TnAUn) ■ --rnXUis) >= L Ct{v) n(ci(A)r (1) 

where Uj is an operator for the primary matter field (tachyon), related to 
top Chern class CT(i^), and r„ is a gravitational operator, related to the first 
Chern class ci of the line bundle Ci at the ith-marked point. We denote 
Tn{Uj) by T„j, and j represents the spin index (j=0,...,p-l). The problem of 
definition (1) has been discussed extensively [11]. 

In a previous paper [18], we have shown that those intersection numbers 
(1) are expressed trhough the correlation functions U (si, s„) as coefficients 
of powers of Sj, 

U{si,S2,...,Sn) = < tre/i^tre/^M . . . ^^^s„M ^ 

II dXie^ < H trS{Xj -M)> (2) 

1=1 1 

where si = iti ; M is sja. N x N Hermitian random matrix. The bracket 
stands for averages with the Gaussian probability measure 

< X >= 1 y dMe->^"+''''^^X{M), (3) 

in which Ais em N xN external Hermitian matrix source. By an appropriate 
tuning of the external source matrix A, we may obtain the desired singularity, 
which generates the p-spin curves. The relation to the generalized Kontsevich 
model is discussed in § 3,4 of [18]. 

An exact and useful integral representation for t/(si, ...,Sn) is known in 
presence of an arbitrary external matrix source A with eigenvalues [19] : 

C/(si,---,s„) = ^(tre^^^---tre^"^) 



N 

N n 



,j + Si 



(4) 
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This representation involves contour integrals around Ui = Oq,. In the 
large N limit, it is convenient to express the factors in the determinant as 
additional integrals. For instance, in the case of the two point correlation 
(n=2), after the shift — >■ — Sj — >■ in the two point function, we 
have 

1 1 

U1-U2+ 2^(Sl + S2) U1-U2- 2^(Sl + S2) 

-/ dxe-<"hh{ — {s, + S2)) (5) 
Si + S2 JO 2i\ 

Tuning now the a^'s, and taking the large N limit, we obtain 

U(si, S2) = ^ -— ^ / dx / duidu2sh{—x{si + S2))e-("^-"^)" 

xexph^ E + ^^^r' - - 2^-0^^^)] (6) 

For the three and four point correlations, similar useful formulae for the 
determinant part of (4) may be found in the appendices A and B of [18]. 



3 Intersection numbers for p = 3 with two 
marked points 

The intersection numbers are obtained as coefficients of the power series in 
Si,S2 of U{si,S2). In a previous paper [18], for p=3, we have computed 
the intersection numbers with two marked points or genus one case (g=l) 
starting from(7). As an example, we compute the p = 3 case up to genus 3. 
The general expansion 

Uis^,S2) = E < rmunrm,,j, >9 r(l - ^)r(l - ^^)sr'^S2^'' (7) 

g,m,j " " 

with the condition, 
(p + l)(2^-2 + n) = ^(pmi + j- + l), m', = m^k + (A; = l,2) (8) 

is applied to the special case n = 2,p = 3 The gamma functions in (7) 
represent the spin factors. 

After rescaling of the parameters, 

Uis.,S2) = ^ J(3..) V3 r ^ysh(^(3.0 V3,) A(. - ^..-/-) 

xM-ay - ^:^,S2'/') (9) 
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in which a = (si/s2)^^^, and the Airy function is 

Ai{y) = J^^ ^e^^+-^ (10) 

The Airy function satisfies the differential equation 

Aliy) = yA{y), A':{-ay) = -a'yM-ay) (11) 

The genus one case (g=l) has been studied in [7]. 

If one expands the hyperbohc sine function and the Airy functions in (9) 
up to relevant orders, we find a sum of six terms which, for = 2, involve 
the following integrals: 

/•CO roo 

h = / dyy^Ai{y)Ai{-ay), h = / dyyA'l{y)Ai{-ay), 
Jo Jo 

roo poo 

h = / dyyAi{y)A'^{-ay), h = / dyyA[{y)A'i{-ay), 
Jo Jo 

h = / dyy^A',{y)Ai{-ay), h = / dyy''A{y)A'^{-ay) (12) 
Jo Jo 

A repeated use of (11) plus integrations by parts allows us to write all these 
integrals in terms of 

3-2/3 I I 3-1/3 I 2 

^•<"> = f(273) = ^'<3)' ^<'°> = -f(T73) = -^^'3) 
plus the integral 

T = / dyA,{y)j^^{-ay) (14) 
Jo 

which cannot be reduced to ^i(O) or ^^(0). For instance one finds 

(1 + a')h = A(0)' - 2T (15) 

and so on. However, all the T-dependence cancels when we sum up all the 
terms relevant to = 2 in C/(si,S2)- For instance the sum of all terms of 

order is given by 



_i^l3i/3(i + a3)2 11 f ^ ^ l.]_^V^a\l + a^sfh (16) 
3! 16 ^ y 2 5 31 V ) 2 b \ ) 

and we add up the six terms and expand in powers of si to the relevant 
orders we find 
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TK., iMO))\ „14/3 2/3 11 11/3 5/3 

U[Si,S2)\g=2 — 32 . 32/3 ^ ^ 2 "5"*! *2 

17 s/3 S/3 11 



< ro,ir4,i >g=2-- 



-y^/^^f (17) 
Prom these results, we obtain the intersection numbers 

1 

11 

< r2,,T2,, >,=2= ^ (18) 

Rather than computing the exact dependence in a of the terms propor- 
tional to S2^^^ and then re-expand in a to obtain the various terms of (17), 
we may proceed in a simpler way by expanding Ai{—ay),A'^{—ay),A'-{—ay) 
for small a: 

A{-ay) = A{0) - ayA'.iO) + jy^A'^O) + ■■■ (19) 

and we then recover (18). 

In the genus- three case (g=3), we have again ten distinct integrals Ji — Jio 
for the terms of order 82(1"^, in the small Si, S2 expansion of (9). 



Jl = 


/ dyy'^Ai{y)Ai{- 
Jo 


-ay), 


J2 = 


roo 

/ dyy'A',iy)Ai- 
Jo 


-ay) 


^3 = 


/ dyy'A,{y)A',{- 
Jo 


-ay), 


J4 = 


roo 

/ dyy'A',{y)A',{- 
Jo 


-ay) 


^5 = 


/ dyy^A'l{y)A,{- 
Jo 


-ay). 




roo 

/ dyy'Aiy)A':i- 
Jo 


-ay) 


J7 = 


/ dyyA':'{y)A{- 
Jo 


-ay). 


Js = 


r dyyA{y)A':\- 
Jo 


-ay) 


J9 = 


roo 

/ dyyA'l{y)A',{- 
Jo 


-ay). 


Jio — 


roo 

/ dyyA',{y)A':{- 
Jo 


-ay) (20) 



The genus 3 contribution for U [31,82) is then expressed as the sum of four 
terms, ^7'-^^ - U^^\ The term U^^\ which is related to Ji, is 



3 



■s^(-15a^ - 42a^ + 90a^° + 63a^^ 



8960 

+630^3 ^ 90^14 _ 42^16 _ i5a^7)Ai{0)A{0). (21) 
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The term U^'^\ which is related to J2, is 
1 



2560 



2560 



sla^^{30 + 72a - UOa^ - 90a^ + 12a^)A(0)^^(0). (22) 



The term U^^\ which is related to J3, is 
1 



C/{3) 



2560 
1 

2560 



sla^{-12 + QOa^ + 1200^ - 720^ - 30a^)^i(0)A^(0). (23) 



The term U^^^ from the sum of the contributions of J4 to Jio- We have 

1 



Jw 
J9 

J7 

J4 



l + a3 
—K2 — Jio 

+ 2a^ - 2a^ - 



(l + a^)^ 



1 



1 + a- 

— nJ& 
a"* 



■J7 



a^ — St/g 



(24) 



with L = Ai{0)A'-{0), and Ki,K2 are given below. We have also the following 
relations between Ji, J 2 and J3, 



J2 
J3 



1 



1 + a- 
1 

l + a3 
1 



:(30J5 + 12J3) 

:-5 J5 + 4 J4) 

;-5aV5 - 4 J4) 



(25) 



Thus U^^^ becomes 



1152 



slia + 5a^ - 20a^ + 23a^° + 16a^^ 



-19a^^ + 130^*^ + 2a^^ + (2a^ + 13a^ - 19a^ 

+16a^^ + 23a^^ - 20a^^ + 50^° + a^^))Ai{0)A[{0). (26) 
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Since a = {si/s2y^^, the above expression for U{si,S2) is a symmetric func- 
tion in si and S2- Denoting 



s"-^"^ = tm,j (27) 

and dividing U{si, S2) by i.e. 1/27 in this case, we obtain the intersec- 
tion numbers < Tm^j^Tmaja >3=3 ^he coefficient of imiji^m2,j2 U{si,S2) 
taking into account the spin factors. The following spin factor appears as a 
over all factor in t/(si, S2) at genus 3. 

■^.0)^5(0) = -(^r(i)r(f ) (28) 

where r(|), r(|) are spin 1 and spin 2 factors, respectively, as (7). All the 
integrals Ji,..., Jio are expressed by (28), and there are no terms like (14), 
which appeared in the integrals for the g=l,g=2 cases. Finally we have to 
compute the following terms 

K^ = J dyA';{y)A{-ay) ^ -A{0)A',{0) - K2 

K2 = [ dyA[{y)A[{-ay). (29) 



For these integrals, we find 

and all the integrals reduces to the spin factor (28). Summing up the results 
of C/^1) to U (4), we obtain the intersection numbers for p = 3, = 3, 

1 1 



< ro,oTv,i >,=3- < ro,ir7,o >,=3- 

_ 5 _ 19 

< Tl,076,l >g=3- ^Yl04' ^ '^1.1'^6,0 > g=Z- 

103 47 

< T2,or5,l >g=3- -7-70C ' ^2,1^,0 > g=2.-- 



217728 ' ' 77760 

443 67 
< r3,or4,i >,=3- < r3,ir4,o >,=3- (31) 

The above results arc in complete agreement with the previous results [13, 
15], which were obtained by recursion relations. 



4 Intersection numbers for p > 3 

For higher multicritical points the algebra is similar, except that we have to 
deal with generalized Airy functions. For instance for p = 4 instead of Ai{x) 
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we have to work with 4>{x) defined as 

POO , . 

(f){x) = / d^;e-4'' (32) 
Jo 

which satifies 

(j)"'{x) = x(j). (33) 

Then, similarly 

3 3 

"l( 1 >l/2,,2 ^2/- 1 U/2,,2 14,14 , , 

where a = (si/s2)^^^- In complete analogy with the p — 3 case, a repeated 
use of integration by parts and of (33) leads to the expansion of U{si, 82)- In 
the genus one case, 

^is,S2f/\s^ + S2){mf (35) 
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with 



0"(o) = 2V2r(^), 0(0) = 2-v^r(i) (36) 

which provide the j = 0, j ~ 2 spin factors respectively. Replacing Si,S2 by 

+ A ('om+(l+i)/p _ J. 

1 3 

U{si,S2)\g=l = 2^^2,0^0,0 + ^l,0^1,0 + ^0,0^2,o)(r(-))^ 

+^(ii,2to,2 + to,2ti,2)(r(^))2 (37) 

Multiplying by a factor ^, we obtain the intersection numbers as coefficients 
of (37) for p = 4 in the genus one case, 

1 1 1 

< To,oT2,0 >g=l= g, < Ti,oTi,o >g=l= -, < To,2Tl,2 >g=l= gg (38) 

18 2 

For g = 2, p = 4, from the term S2 Si , we have similarly 

< ro,ir4,i >g=2= ^ (39) 

For general p we have to deal with the generalized Airy functions (f){x) 
for p > 2, which satisfy the differential equation, 

(f)^\x) ^ x(f){x) (40) 
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where (f)^^\x) means the p-th derivative of (j){x). The generahzed Airy func- 
tion has an integral representation, 



(f)(y) = / du e--+^". (41) 
Jo 

As examples of what the method can provide we give a few results : for the 
case p = 5, we obtain 

1 

< T-i,3ro,2 >9=1 = < Ti,2ro,3 >s=l= ^ 

1 

< ''■l,0'?'l,0 >g=l = < To,oT2,0 >fl=l= g 

< ro,,n, >,.2= (42) 

For the case p — 6, 

1 1 1 , , 

< ^0,3^1,3 >g=l= ^, < T-o,2ri,4 >g=l= — , < ro,4ri,2 >g=l= — • (43) 

For the case p — 7, 

1 

< ro,2ri,5 >g=l = < Ti,2To,5 > g=l= ^ 

1 

< To,4Ti,3 >3=1 = < To,3Ti,4 >p=l= — 

< Ti,ori,o >g=i=< ro,or2,o >s=i= J- (44) 

5 The dependence of the intersection num- 
bers 

In a previous article [18], we have considered the intersection numbers with 
one marked point for arbitrary p, and found results such as 

p-1 



< Tin >g=l 



24 

P' 



_ (p-i)(p-3)(i + 2p) r(i-^) 

'=^~ p • 5! • 42 . 3 r(l - i±i 



p ' 

{p - 5)(p - 1)(1 + 2p)(8p2 - I3p - 13) r(l - 



< ^^=3 p2 . 7! . 43 . 32 _ i±i) 

(p -7)(p- 1)(1 + 2p)(72p4 - 298p3 - ITp^ + 562p + 281) 



< T^ J >g=4 



p3 . 9! . 44 . 15 

r(i--) 
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with n = 2g — 1 + ^ ^ . In the large p limit, the intersection numbers 
< Tn,j >g behave as 

with Bg is a Bernouilli number, Bi — 1,^2 = -^3 ~ -^4 ~ Note the 
well known relation to CC^q) ^ 

(Mm - (^«^*' 

We have derived (46) from U{s) in the large p limit. The one point 
function U{s) has the following expression [16], 

With s = -, and = x^, we have 

Thus we obtain 

This may be written as 

- l'^l^^^r<rhM-ho,f) (51) 

and expanding the exponent in i, we find 
If we use the expansion 

sh- °° R n-^" 

and drop the factors (^^)~^, and {a/2)~^^^ which are close to one in the 
large p limit, we obtain 

^(s) = (1 - - E(-l)"~'77r^^'")i^r(l + -) (54) 
^ ^ ^ p^i^ ^ (2n)!2n ^ Na ^ p' ^ ' 
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Since the intersection numbers < r^j >g are related to U (s) by [16] 

U{s) = E < rn, >, ^r(l - i±l).(^-^)(^+ly- (55) 

with {p + l){2g — 1) = pn + j + 1, we have rederived the large p behavior of 
(46). 

Prom (53), taking a derivative with respect to cr, gives, 

111°° R 

;;^4-^i:(-ir-|^^-'. (56) 

Using this relation one obtains 

±(,U(a)) = \-\-^^ (57) 
The di-gamma function ■0(^) has the following expression, 
H.) = A,ogr(.) = log. - 1 - f <i.(i - i + -L.)e^". (58) 



From (57) and (58) we find thus in the large p limit, 

A,ogr(.)= ,„g,_^ + ^-,,(^(.t/(.)))e-" (59) 



= log2; z — / daU(a)e 

2z dz J 



The last integral is related to the density of states. In (2), s is replaced by 
s = it, and if we replace z by iE, and take the imaginary part, we obtain the 
density of states p{E). After integration by parts, we obtain 

p(£;) = ^imiogr(i£;)-^-^ (6o) 

We will discuss this expression in connection to the density of states of 
the SL{2, R)/U{1) black hole sigma model in the next section. 

Next we consider the two point correlation function U (si, 52)- For general 
p, C/(si,S2) is expressed as 

^(^1,^2)= . , w TYT- / dx dvidv2sh{ — - — (psi)^/^) 

[Si + S2)[pS2)'^'P Jo Jo 2 
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The exponent of (61) follows from the binomial expansion, 

(u + -^r^' = u^-^' +ip+ i)u^C-) + ip+ i)p\u^-\'^r +... (62) 

and we use c{p + 1) = 1, vFs = v^/p. As in the case of p=3, polynomials 
in a (21) give the intersection numbers. Therefore we expand (61) in power 
series of a. At lowest order in a, we obtain two terms from (61), 

3 

C/i = 7^(51 + 82)'^^^ I dxx^ct>{x)ct>{-ax) 

U2 = - ^^^~^\ -)'sl{ps2)'-^ [ dxx(t>{x)4>^-''\-ax){-af-^ (63) 
24 S9 J 



2+i- 

Prom U2 we find a term proportional to 0^2 ^ , namely 



AU2 = ^pia4^'(</.(^-2)(0))^ (64) 



with 



Jo 



= p"pr(i--). (65) 
p 

2+i i 

Since S2 Si = ^2,0^0,0; we obtain 

< To,oT2,0 >g=l- (66) 

2+- 

Prom C/i and U2, we collect terms proportional to a^S2 " and obtain 

< To,2Tl,p-2 >s=l= (67) 

This result agrees with those obtained previously for p = 4, 5, 6 and 7 in 
(38), (42), (43)and (44). The intersection number (67) can be neglected in 
the large p limit in comparison with (66) . 

Similarly one obtains the g=2 terms from the coefficients of a"^S2 ^ 
(m=l,2,3), 

(p-l)(p-3)(2p+l) 



< To,oT4,2 >ff=2 = 

< 7o,l7"4,l >g=2 — 



5760p 
(p-l)(p-2)(p + 2) 
2880p 



. ^ (p-l)(p-3)(2p+ll) 

< ro,2r4,o >.^2= (68) 
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For the particular values of p = 3, 4, 5, the above expressions agree with the 
previous results (18), (39) and (42) for the genus two case. 

4+4 

From the a^S2 ^ term, one finds 

2p3 + i3p2_i58p + 215 
< ro,4r3„-2 >,.2= (69) 

which is valid for p > 6. 

In the large p limit, the three terms of (68) become equal, and coincide 
with the result for the one point intersection number (46). 

< To,mT4,2-m > g=2= (p ^ Oo) (70) 

Note that (69) is order p, and is negligible compared to (68). 

From the terms a™S2 ^ in the small a expansion of U{si^S2)i we obtain 
the g=3 (genus 3) terms. In the case m=l, we have 

^ ^ (p-l)(p-5)(2p+l)(8/-13p-13) 

This is identical to < r^^^ >g=3 in (45). The identity follows from the string 
equation, in which the insertion of ro,o reduces the intersection number from 
s to s — 1 marked points: 

s s s 

< To,0 n Vji >3= Yl < ^ni-l,ji n ^nuji >g (72) 

i=l 1=1 i=l,i^l 

In our formulation, this string equation follows from the integral represen- 
tation for the intersection numbers, when one collects the terms proportional 
to a. By explicit calculation of two marked points, we verified this string 
equation. It might be possible to verify this string equation for n-marked 
points by the taking account of the term of a. 

Prom a S2 " , we have for p > 5, 

^ ^ _ (p-l)(p-2)(p-4)(p + 2)(2p+l) 

< To,lT6,3 >g=3- ^2 . 7! . g . 32 y''^) 



From a^s 



_ (p - l)(p - 3)(16j/ + 34p^ - 155p - 129) 
< 7-o,2r6,2 >5=3- p2 . 7! . 64 . 32 ^'^1 

In the large p limit, these g=3 terms exhibit same behavior as in (46), 
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6 Analytic continuation to negative p 



One may analytically continue the integral representations of the correlation 
functions to negative values of p. This continuation was already examined 
in [18], and we recall some of the results here : 

U(s) = — / i!ie-^[("+^^)^^'-("-5^)''^^l (76) 
^ ' Ns J 2m ^ ' 

where c= ^E^- 

Expanding the exponent, we obtain 

= / ^exp[-c(.«>+j^.3„,-.^ P(P - 1)(P - 2)(P - 3) ^.„,_.^, . .J, 

(77) 

This integrals yield Gamma functions after the replacement u — (— )^/^, 



^ ^ Nsp {csy/p Jo 

1 1 rP-l^rn ^ , (p-l)(p-3)(l + 2p) 3 

cpi/r(i--) + ^, ,n o i/r(i--) 



A^sp (cs)Vp^ 24 ^ ^ 5! -42 -3 " ^ p 

{p - 5)(p - 1)(1 + 2p)(8p^ - 13p - 13) 3 _ 5 
714332 ^ ^ p' 

+ (P - 7)(p - 1)(1 + 2p)(72p^ - 298p2 - 17^^ + 562p + 281) 

with y = cps'^'^p. 

From this expansion, we obtain the intersection numbers for one marked 
point as (45). The intersection number < Tnj >g is obtained from the term 
y^r(l-i-i)in(78). 

The continuation to p < is straightforward. The integral in (78) can be 
changed to u by t = ^ ,(0 < f < 00), and one obtains the small s expansion for 
negative p. Therefore the expression for the intersection numbers (45) can be 
analytically continued to negative p. This analytic continuation can also be 
done for two marked points, since we have computed them in the previous 
sections for general p. For instance, from (45), we have the intersection 
numbers for p — —3, 

_ 1 _ 1 

< Ti,o >fl=l- -g, < T3,2 >g=2- ^ 

< re,i >.^3= (79) 
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In a previous article [18], we have computed the intersection numbers 
< Ti,o >g for the case of p = —1 from U (s), which provides the orbifold Euler 
characteristics x(A4g,i) with one marked point, 

< ri,o x{Mg,i) = C(l - 29) = (80) 

with the Bernoulh number Bg, {Bi = ^,-82 = ^^-^s = ^,---)- The s- 
point orbifold Euler characteristics x{-^g,s) may be obtained from the dilaton 
equation: 

< n,or„i ji ■ • ■ r„, J, >g= {2g-2 + k)< ■ ■ ■ r„, >g (81) 

Since the Euler characteristics with s marked points is < Ti^q " " 'Ti,o >g, the 
dilaton equation yields from (80), 

x{Mg,s) {n,oy -^1^(2^ + « - 3)!^. (82) 

This agrees with previous results obtained in [20, 21, 22]. 

For p — —2, we have considered previously the equivalence with the 
unitary matrix model in a matrix source [26] . 

The central charge of the gauged Wess-Zumino-Witten model with sym- 
metry SU{2)k/U{l) is 

^ = '-kT-2 ^''^ 

Changing pto p— — p', k to k — —k' {p < 0,k < 0), we have p' — k' — 2, and 
the central charge C is given by 

C = 2 + jA_ («4) 

The analytic continuation to negative p yields a gauged WZW model for 
S L{2 , R) k' / U (1) . It is known that this model represents a black hole a 
model[9], in particular for the value /c' = | (p = — |), for which the central 
charge C becomes 26. 

The density of states for the SL{2, R)/U{1) black hole has been studied 
in [24, 25, 27], 



1 , 1 , r(-/^ + i - //Or(-/^ + - + m. 



^^^^ = n'""^' + A^^dE'''^ ^{+^E+l + rhn+^E+l-rn) ^''^ 

in which e is a regularization factor, and m — |(n — kw), fh = —^{kw + n) 
are eigenvalues of Jq and Jq in CFT (J^ - ^ n, Jq + Jq ^ -kw). If we 
neglect m , m, and the | terms in the large E limit, we obtain 

1, 1 , r(-iE) 
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or 



p{E) 



2 d 
ndE 



Im logr(-i^) 



(87) 



This expression agrees with (60), obtained from the intersection numbers 
for large p. We have scaled s = cr/p, and the expression (60) is valid for 
small s. Therefore, the Fourier transform of U (s) gives the large E behavior, 
in which the terms m,m and 1/2 in (85) can be neglected. 



In this article, we have shown that the correlation functions U{si, S2, - ■ ■ , Sn) 
of a Gaussian matrix model in a tuned external source, provide the intersec- 
tion numbers for p-spin curves. For instance, from the two point function 
U{si,S2), in the case of p=3, the intersection numbers are computed up to 
genus 3, 

We have also computed the intersection numbers for general p. They are 
given by power series in a, a = (f^)^- Then we have considered the large p 
behavior for the two point functions. The density of states p{E) becomes a 
di- gamma function in the large p limit, and this expression agrees with the 
density of states of a SL{2, R)k/U{l) WZW model, which has been studied in 
the context of two dimensional black hole solutions. The n-point correlation 
functions U{si, ■ ■ ■ , s„) are known through the determinant of a kernel for the 
p-spin curve case. It will be interesting to investigate further the detailed 
comparison of those correlation functions, between SL{2, R)k/U{l) WZW 
theory and the intersection numbers for negative p-spin curves . 
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